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Constraining dark energy with gamma-ray bursts 

Lado Samushia^'^ and Bharat Ratra^ 

ABSTRACT 



We use the measurement of gamma-ray burst (GRB) distances to constrain 
dark energy cosmological model parameters. We employ two methods for analyz- 
ing GRB data — fitting luminosity relation of GRBs in each cosmology and using 
distance measures computed from binned GRB data. Current GRB data alone 
cannot tightly constrain cosmological parameters and allow for a wide range of 
dark energy models. 



Subject headings: cosmological parameters 
ture of universe — GRB 



distance scale — large-scale struc- 



1. Introduction 



The combi nation of recent measurernents of distant supernovae Type la (SNe 



ent mag nitudes JSahni et al.lboOsilCunhJboogl J iPerivolaropoulos fc ShafielodI 



20091 ) , cosmic microwave background anisotropy flDunklev et al 



the baryon acoustic signa. 



2005 




Percival et a 


(Allen et al. 


2008: 



2007 



in the power spectrum of ga laxies ( lEisenstein et al. 



Samushia &: Ratra 



2008 



200 



2009 



a) appar- 



Hicken et al 



Komatsu et al 



2005 



20091 



Cole et al. 



Samushia &: Ratrall2009al). and gala xy cluster gas mass fractions 



Ettori et al.ll2009l ) indicates at high confidence 



that about 70% of the energy budget of the universe comes from non-luminous dark energy 
that is close to spatially uniform and has negative pressure 
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therein. For a recent review see 



(2008) 



Silvestri fc TroddenI (|2009l ) 



Bamba et al. ( 2009 ), Aluri et al] ( 2009 ). and references 
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Although the existence of dark energy is now a well-established observational fact its 
physical nature is still a topic of great debate. The simplest and historically first physical 
model of dark energy is the cosmological const ant A that ha. s an equation of state p = —p, 
where p is its pressure and p the energy density ( lPeebleslll984( ). The cosmological constant is 
introduced by hand into the equations of general relativity but could be related to the energy 
of the fluctuating vacuum. This ACDM model, although sim ple and in good accord with 



most available cos mological data§ (e . g., 



Frieman et al. 



20081 ). has a number of theoretical 



shortcomings (e.g.. iRatra fc Vogeleyl l2008l ). These include the so-called "smallness" and 
"coincidence" problems. The measured value of the cosmological constant energy scale is 
of order 10~^eV, some 30 orders of magnitude smaller than the Planck scale, perhaps what 
we would expect the vacuum energy scale to be based on simple quantum field theoretical 
arguments. In addition, since the energy density of nonrelativistic matter is decreasing with 
the expansion of the Universe while the energy density of a cosmological constant is constant, 
there is a very narrow window in time when both of these have comparable energy densities, 
and it is unclear why we happen to live at this special time. 

Because of these and other issues a number of different dark energy models have been 
considered. Typically these introduce a new component that acts like a cosmological con- 
stant, in that it is close to spatially homogeneous, while gradually decreasing in time to the 
small currently observed value. An early example of this class of models is the (/)CDM model 
in wh ich dark energy is taken to be a scalar field ( iPeebles fc Ratralll988l : iRatra &: Peebles 
19881 ). In this model, the current cosmological constant energy scale is small because th e 
Universe is old . For recent revi ews of dark energy see ICaldwell &: Kamionkowskil ( 120091 ). 
FriemanI toO^ . andlSami J^OoJE 



Although available cosmological data can constrain dark energy models, they are not 
yet good enough to strongly discrimiriate between dif ferent dark energy models (see, e.g.. 
Gong et al.l l2009l iKilbinger et al.l l2009l ICoc et all l2009l and references therein) . In the near 
future better-quality and more ii idependent cosmo l ogical observations should be able to 



break this degeneracy (see, e.g., IWang et al.l l2009l : iMortonson et al.l l2009l : iThomas et al. 



^The ACDM model assumes the cold dark matter (CDM) model of structure formation which might have 
some observational inconsistencies (see, e.g., Peebles fc Ratra 2003t Perivolaropoulos 20091 : Primack 20091 
and references therein). 

•^In the 0CDM model, we consider in this paper the dark energy scal ar field couples t o the matter only 



through gravity. For da r k ener g y mod e ls wit h less restrictive c oupli n gs see I Wang fc ZhansJ (|2008l ) , iDent et al 



(|2009f ). lLa Vacca et al 



(200 



i), 



J ami 



Chongchitnan J 2009 ^ Nesseris J 2009 ), a nd ref er ences therein 



For o t her dark energy niodels see [Bihc et al. (l2009h. iBasilakos fc Perivolaropoulos! ( 20081 ). Grande et al 
I2OO9I ) , butta k Scherreij |2009l) , li^ng (2009, ) , Undrianov et al.l (|2010l ). and references therein. 
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20091 : lArun et al.lbood : IVashar et allboOoh FI 



SNe la were one of the first cosmological probes to give direct evidence for dark energy. 
These are very bright exploding stars and are standardizable candles that can be seen to very 
large distances. More then 300 well-calibrated distant SNe la h ave already been observed up 
to the redshift of 1.6 (IKowalski et al.ll2008l : iHicken et al.l 120091 ). These measurements alone 
give more then Sex evidence for the existence of dark energy, but they are not very effective 
in constraining dark energy model parameters overall unless used in combination with other 
data. Even when current SNe la data are combined with all other currently available data 
sets, we cannot yet determine if the energy density of dark energy is constant as required 
by the ACDM model or if it varies in time as suggested by dynamical alternatives such as 
(pCDM. A dedicated SNe la space mission shoul d result in significantly more higher qu ality 
data that should help resolve this issue (see, e.g.. iPodariu et al.ll2001al : lAlam et alll2009l and 
references therein.). 

One way of improving our understanding of how dark energy behaves is to study the 
evolution of the universe at redshifts higher than those probed by SNe la. This requires 
standard candles that are visible at great distances. Gamma-ray bursts (GRBs) could in 
principle serve as such high redshift standardizable candles. They ar e the most lumino us 
events in the universe today and can be seen beyond z = 8 (see, e.g., iTanvir et al]l2009h FI 
If it is definitely established that GRBs are standardizable candles, their visibility at high- 
redshift should prove to be very useful in discriminating between ACDM and time-varying 
dark energy modelsj^ 

With the intention of getting cosmological const raints from GR B observations a number 
of GRB calibrations have been used so far (see, e.g.. ISchaefed 120071 . and references therein). 



In the near future, measure ments of nonlinear structu r e formation wil l help d i scriminate between diff erent 



dark energy models, see, e.g., iGrossi fc Springel ( 20091) . Francis et al. ( 2008 ). Casarini et al. ( 20091) . and 
references therein. Other tests that h old significant p o tential are the angul ar size of radio sources and 
galaxies as a function of redshift (e.g., Daly et al. 2009 : Santos &: Lima 20081 and references therein) and 



2008; 


Dev et al. 


20081 


Fernandez-Martinez & Verdel 


2008 



^For a review of GRB physics see, e.g., Meszaro^ (2006). 

^For early d i scussions of the use o f GRBs as a cosmology pr o be see, e.g., Lamb fc Reichar t' (200(f), 



Nemirofi (2000). 



Ghirlanda et al 



(2004). Friedman fc Bloom (2005), 



Firmani et al 



Xu et al..(.2005.) . 



Mortsell fc SoUerma n (20051). lOi Girolamo et al.l J2005l). iBertolami fc Silval 12006^ a n dlLamb et al.l ( 2005 ). 



More recent studies may be traced back through iMosauera Cuesta et al.l (120081) . lAmati et al.l (l2008l) . 



Basilakos fc PerivolaropoulosI (|2008l ). ICapozziello fc Izzol |2008l) . iTsutsui et all |2009l) . and lOi et all (|2008l ). 
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One that gives least scatter and therefore most information is 



log 



1 erg 



Ai + Bi log 



peak 



l + z) 



300 keV 



a relation that connects the total burst energy of the GRB {E^) to the peak energy of 
the GRB spectrum (-Epeak) (IGhirlanda et al.l l2004j ). Regrettably, we do not yet have a 
model-independent way of computing the coefficients Ai and Bi. A better understanding 
of physical processes that result in the burst, or observations of nearby GRBs (to which 
distances can be measured independently), could in principle help us to calibrate the E^- 
Epeak relation without any prior assumptions. To extract cosmological information, GRBs 
have to be recalibrated for every dark energy model considered (at each set of parameter 
values). This is time consuming and also results in large statistical uncertainties and hence 
GRB cosmological constraints that are poor. 

Recently, methods of calibrating GRBs in cosmology-independen t manners have been 



proposed and used to constrain sorne dark energy models (see e.g., iKodama et al.l 12008 



Liang et al.ll2008t IWei fc ZhangI l2008l : iLiang fc Zhangll2008l . that use SNe la measurements 
to externally calibrated GRBs). The resulting cosmological constraints are still loose, but in 
the future when more high precision GRB observations become available this could provide 
a strong test of dark energy. 



Wangl ( 120081 ) recently used data of 69 GRBs (ISchaeferl 120071 ) to construct a distance 
measure that can be used to constrain cosmological models. The advantage of this method 
is that internally calibrated GRB data may be straightforwardly combined with other data 
when deriving cosmological constraints. On the other hand, with this method the resulting 
cosmological results are sensitive to the chosen binning. This method also requires an input 
cosmological model and thus is not completely cosmology independent. When this method 
is used to constrain ACDM the GRB data favor lower values of both cosmological constant 
energy density (^a) and nonrelativistic matter energy density (fim) than do the SNe la data. 
The GRB data by themselves are unable to strongly constrain cosmological parameters, for 
example in spatially flat ACDM the GRB data require flm 
JWaneJl2008h . 



0.25^0;^^ at la confidence 



In this paper we use GRB data to constrain time variation of dark energy's energy 
density. First we recalibrate GRBs for each cosr nological mod el and compare the result with 
the ones derived using the data and method of IWangj (120081 ). We consider a 0CDM model 
where a scalar field which is close to spatially uniform on cosmological scales slowly rolls 
down an almost fiat potential and plays the role of dark energy. 



In the next section, we summarize the dynamics of this scalar field dark energy model. 
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In Section 3, we describe our methods and computations. We present and discuss our results 
in Section 4. 



2. Scalar field as dark energy 

In the 0CDM model, consi stent with the indicat i ons from cosmic m icrowave background 
anisotropy measurements (e.g.. iPodariu et al.ll2001bt iPage et al.ll2003l ). we only consider the 
spatially flat universe case. The invariant four-interval in a homogeneous and isotropic 
version of such a universe is 

ds"^ = -dt^ + a?{t)dx ■ dx, (2) 

where t is cosmic time, dx is the spatial separation in three-dimensional Euclidean space and 
a{t) is the time-dependent scale factor which determines the change in distance between two 
distant noninteracting test particles in the universe. 



In the set of models we consider the scalar field with Lagrangian density 



Id, 



(3) 



is the dark energy. Here V{(j)) is the potential energy density of the scalar field. In the 
expanding, spatially homogeneous and isotropic universe described in Equation the spa- 
tially homogeneous scalar field obeys the modified Klein-Gordon equation. 



+ 3-0+- = 0, 

a 2 O0 



(4) 



and the dynamics of the scale factor is governed by 



SttG 



(Pm + 



(5) 



In Equation ([5]), G is the universal gravitational constant, and pm and p,^ are the energy 
densities of nonrelativistic matter and the scalar field, respectively. If the scalar field is 
uniform in space then, from Equation ([3]), the energy density of the scalar field is 



P<^ 



327rG' 



0' + V^(0) • 



(6) 



Since an underlying more fundamental explanation of dark energy remains elusive, there 
is as yet no first principles way of choosing the scalar field potential. If we choose the 
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potential energy density to be inversely proportional to a power of t he scalar field, V((t>) 



, the (t>CDM mode l has a number of very interesting features ( jPeebles fc Ratral Il988 



Ratra &: Peebles! Il988l ). First, even if the scalar field starts off from a very high energy 
density state, its energy density decreases to a very small value during the course of cosmic 
evolution. Second, in the radiation and matter-dominated epochs the evolution of the scalar 
field "tracks" the evolution of the dominant component. The scalar field slowly comes to 
dominate, leading to the end of matter domination and the start of the scalar field-dominated 
epoch. So, in the 0CDM scenario the "smallness" and "coincidence" problems mentioned 
above are partially resolved because of the time-evolution properties of the scalar field. Even 
if the scalar field potential is not exactly inverse power law, this form of potential provides 
a very economic way of parameterizing the slowly-evolving dark energy scenario with just 
one positive parameter a. Moreover, unlike the XCDM para meterization of dark energy, the 
0CDM model is physically consistent (see, e.g., lRatralll99ll ). 



In the 0CDM model, the dark energy density, unlike the cosmological constant, varies 
slowly in time. The dark energy density increases as we go back in time and larger values of 
a correspond to faster evolution of dark energy. As a result, observable quantities such as 
luminosity and angular diameter distances in the 0CDM model differ from the predictions 
of the ACDM model with a time-independent cosmological constant. Some predictions, of 
course, depend also on the values of other cosmological parameters, such as the Hubble con- 
stant Hq or the energy density of baryonic matter fib? and so many independent cosmological 
tests are required to break this degeneracy and constrain dark energy model parameters. 

A number of cosmological tests have been used to constrain the 0CDM model, includ 



ing the angular sizes of radio sources and q uasars as a function of redshift ( IChen fc Ratra 



2003: Podariu et al. 



tion of redshift ( Wilson et al 



2003; 



Da: 



y et ahl 120091). the apparent mag nitude of SNe la as a func 



20061 : ISamushia fc Ratral l2009br). and the gas mass fraction 



of lar ge relaxed clusters as a function of redshift (jChen fc Ratral l2004l : ISamushia fc Ratra 
20081 ) . Current cosmological observations are in good agreement with a time- independent 



cosmological constant in a close to spatially flat ACDM model, but slowly rolling scalar fleld 
dark energy in the 0CDM model is not yet ruled out at high confldence. 

In this paper we use measured GRB luminosity distance as a function of redshift to 
constrain slowly rolling scalar fleld and other dark energy models. Because of the calibra- 
tion problems mentioned above, GRB data alone cannot constrain cosmological parameters 
effectively. We also combine the results obtained from the GRB analysis with constraints 
from SNe la and baryon acoustic peak measurements, to illustrate the effect and weight of 
current GRB data in such a combined analysis. 
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3. Cosmological constraints from GRB 



Besides ttie Epeak^E^ relation, (Equation ([T])), ISchaeferl (120071 ) uses four other cali- 
brations for GRBs that relate total luminosity to directly measurable quantities. These 
calibration relations are given by 



log 
log 
log 
log 



lergs^^ 
L 

lergs^^ 
L 

lergs^^ 
L 

lergs^^ 



A2 + B2 log 
A3 + 53 log 
+ Bi log 
+ 55 log 



^lag(l + Z 

0.1 s 

V{l + z\ 



1-1 



0.02 



E. 



peak 



1 + ^) 



300KeV 



(7) 



0.01s 



where L is the absolute luminosity, rjag is the GRB lag time (the time shift between the hard 
and soft curves), V is the variability (the normalized variance of an observed light curve 
around the smoothed light curve), -Epeak is the peak energy of the GRB, and trt is the rise 
time or the time over which the light curve rises by half of the peak flux. 



Following ISchaeferl ( 120071 ) we take Equations ([7]) and ([T]) and for each cosmological 
model find the best-fit value s for the A and B parameters using the bisector least-square 
method (for a description see llsobe et al.lll990l ). We use the best-fit A and B values and the 
measured riag, V , -Epeak, trt values to compute L and E^ using the same Equations (j7]) and 
([1]). For each calibration relation we then compute the luminosity distance as 



dl 



dl 



L 



47rPboio 
^-^(1 + z) 

4vrS'bolo-^beam 



where Pboio is the bolometric peak fiux and Sboio is the bolometric fiuence of the GRB. Fbeam 
is so-called beam factor Fbeam = 1 — cos6'jot, where ^jet is a jet opening angle. 

We then derive an effective luminosity distance by weighting the five estimates of GRB 
luminosity distances. 



log(<(^.)) 



Ealog(4(^^)«)/<a 



a^(log(4(^.))) = l/EV^. 



2 



(9) 
(10) 
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where index i runs over 69 redshift bins, and a over five calibration relations. To constrain 
cosmological parameters we use defined by 



2 _ (log(4(^.))°^^ - ^ogidli^^)rr 



We also adopt the method of lWand ( l2008l ) for using GRB data to constrain cosmological 
parameters. S he placed each of the 69 GRBs in the redshift range z = 0.17 to ^ = 6.6 
( ISchaefen 120071 ) at a luminosity distance that minimized a combined that took weighted 
account of all five calibration relations. She then computed a distance measure 



rp(0.17)' 



where 



Ho 



(12) 



(13) 



/IC2(1 + Z)V2' 

and dL{z) is the luminosity distance at redshift z, h = Hq/{100 km s~^ Mpc~^), and c is the 
speed of light. The ratio in Equation f lT^ does not depend on the Hubble constant and does 
not require information about the absolute calibration of GRBs (which are unknown). 



Wand ( I2OO8I ) computed the distance measure in six redshift bins 



rp(zi), i = 1,2,...6. 

The values of fnfzj) are shown in Table II and the normalized covariance matrix is shown 



in Table III of [Wang (120081 ). For currently viable cosmological models, these rp{zj) are 



almost completely inde pendent of t he cosmological model and so provide a useful summary 
of current GRB data (jWangj |2008| ). This information can be used to constrain any dark 
energy model and the resulting GRB data constraints can be straightforwardly combined 
with other constraints. In this approach, is given by 



x\nra,p) = A{zO(r,{S- 



where 



A(zO 



(14) 



(15) 



Sij is the normalized covariance matrix given in Table III of [Wang (120081 ) and summation 
over repeated indexes is assumed. Here, cXi is if A{zi) > and if A{zi) < 0. 

In this paper we consider three cosmological models, ACDM, the XCDM parameteri- 
zation of the dark energy equation of state = oj^Px in a spatially flat universe, and the 
spatially fiat (pCDM model. Since we are comparing low-redshift predictions to observations 
we ignore the contribution of radiation in these models. In this case, in all three models, 
the background evolution can be fully described by two parameters, the fractional energy 
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density of nonrelativistic matter and a parameter p that describes the properties of dark 
energy. In ACDM p is the fractional energy density of the cosmological constant Q\, in 
XCDM it is the equation of state parameter u^, and in 0CDM it is the positive parameter 
a that governs the steepness of the scalar field potential energy density. 

In each model, for both methods, we divide the two-dimensional space of cosmological 
parameters into an equidistant grid, and for each pair of parameters flrn and p we compute 
the theoretical luminosity distance. We then compute the difference between the t heoretical 
predic tion and the measured value at each of the 69 redshifts listed in Table 2 of ISchaefer 
feoOTh or the six redshifts listed in Table II of lWaneJ fcoosh . 



The best-fit parameters are defined as the pair that gives the minimum value 

of x^{^m,p)- The la, 2a, and 3a confidence level contours are defined as the sets of points 
where the value of x^i^m^p) is more than its minimum value x^i^m^P*) by 2.30, 6.18, and 
11.83, respectively. If the likelihood function, oc exp (— x^/2), was Gaussian, the 3a contour 
would be 99.7% likely to enclose the true values of cosmological parameters. In our case 
the likelihood function has a single maximum and decreases monotonically from the best-fit 
value point, so the true values of cosmological parameters are very unlikely to be outside the 
3a contours we compute. 



Figures 1-6 show the constraints from the IWand fl2008f ) GRB data on ACDM, XCDM, 
and 0CDM mod el parame t ers. F igures 7-12 show the constraints derived using a method 
similar to that of ISchaeferl (120071 ). 



Results and discussion 



The GRB constraints shown in Figures 1-3 are consistent with the "standard" spatially 
fiat Qjy = 0.7 ACDM model at a little under 2a, with the GRBs mildly favoring a somewhat 
lower value of than the "standard" value that is compensated by space curvature and 

= (Figure 1) or a mildly time- v arying dark energy (Figures 2 and 3). While the 
constraints of Figs. 1-3 derived from the IWangl (120081) data ha ve not previously been shown, 
these results are impl icit in the disc u ssions of IWana (120081) and o t her analyses o f GRB 
con straints, including iKodama et al.l (120081 ). iLiang fc Zhangl (120081 ). lAmati et al.l (120081 ) . 
and iTsutsui et al.l (120091 ) . 



These results, and the fact that the current GRB data contours are quite broad, are prob- 
ably an indication of the preliminary nature of the current GRB data. Current GRB data 
by themselves are unable to effectively constrain cosmological parameters. The GRB con- 
straints however, are a little tighter than radio galaxy angular size versus redshift constraints 
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(e.g., iDaly et al.ll2009l ). constraints from strong gravitational lensing data (e.g., iChae et al. 
2004 ) , and those from the measurement of the Hubble parameter as a function of redshift 



[e.g., 



Samushia et al.l 120071 ) 



To get tighter constraints on cosmological parameters, and to see how current GRB data 
affect constraints derived from other data sets, we combine the results of ou r GRB analysis 



with SNe la apparent magnitude versus redshift data (the Ur iion data of iKowalski et al. 



20081 ) and measurements of the baryon acoustic (BAO) peak (jPercival et al.l 120071 ). Since 
all three sets of measurements are independent, we define the total likelihood function of 
cosmological parameters for the combined data as a product of the individual likelihood 
functions. 



Aot — '^^GRB'^SN'^BAO! (16) 

and compute the best-fit values and confidence level contours from £tot as before. The SNe 
la likelihood function Csn depends on the assumptions that we make about the value of the 
Hubble constant. Here we marginalize over h = 0.73 ± 0.03 with the 

4-6 show the constraints on cosmological parameters of the ACDM, XCDM, and (pCDM 
models from a joint analysis of the SNe la Union and baryon acoustic peak measurements, 
without and with the GRB data. These plots show that current GRB data only marginally 
affect the joint SNela and BAO peak constraints (which are amongst the tightest provided by 
current data), favoring slightly lower values of the nonrelativistic matter density parameter 
^Ill- 
Figures 7-9 are similar to Figures 1-3, b ut derived by r ecalibrating GRB data for each 



cosmology, using a method similar to that of ISchaeferl (120071 ). The "standard" spatially fiat 
= 0.7 ACDM model is about 2.5a from the best-fit value. In all three models GRB data 
favor a nonrelativistic matter-dominated universe. Figures 10-12 show the joint constraints 
from GRB, SNe la and BAO data. The constraints are dominated by the SNe la and BAO 
data, with the GRB data shifting the best-fit values to slightly larger values of fim- 

The GRB constraints derived using the two different methods disagree with each other 
at more than 2a confidence level (compare Figures 1-3 and 7-9). This is somewhat worrying 
but not completely unexpected as the field is still under rapid development. 

GRB data alone do not provide tight constraints on cosmological parameters. Moreover, 
while not greatly significant, current GRB data favor cosmological parameter values that are 
at odds with what other data favor. When used in combination with some of the highest- 
quality current data (e.g., SNe la and BAO peak measurements) current GRB data only 
slightly change the results. This is mainly because in the absence of an absolute calibration 
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of GRBs they, as standard candles, have big measurement uncertainties. This is however 
quite hkely to change as more and better-quahty GRB measurements become available with 
improvements in methods to calibrate GRBs. GRBs could potentially provide a very strong 
test of the time variation of dark energy as they can be observed up to redshifts beyond 
eight, at distances where other standard candles cannot be detected. 

Wc thank the referee for a detailed and helpful report. We acknowledge support 
from DOE grant DE-FG03-99EP41093 and the Georgian National Science Foundation grant 
ST08/4-442. 
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Fig. 1. — la, 2(7, and 3a confi dence level co ntours for the ACDM model from the GRB data, 
derived using the method of IWangI (120081 ). The circle indicates best-fit parameter values 
fim = 0.16, Qa = 0.0 with = 0.41 for 4 degrees of freedom. The dashed line demarcates 
spatially flat ACDM models. 
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Fig. 2. — la, 2a, and 3a confi dence l e vel co ntours for the XCDM model from the GRB data, 



derived using the method of IWangI (120081 ). The circle indicates best-fit parameter values 
Qra = 0.0, = —0.52 with = 2.17 for 4 degrees of freedom. The dashed line demarcates 
spatially fiat ACDM models. 
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Fig. 3. — 1(7, 2(7, and 3a con fidence level contours for the 0CDM model from the GRB data, 
derived using the method of IWangI (120081 ). Numerical noise is responsible for the jaggedness 
of some parts of contours. The circle indicates best-fit parameter values Q^a = 0.0, a = 10.2 
with = 1-39 for 4 degrees of freedom. The a = horizontal axis corresponds to the 
spatially flat ACDM case. 
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Fig. 4. — la, 2a, and 3a confidence level contours for the ACDM model (dashed line demar- 
cates spatially flat models) . Solid lines (circle denotes the best-fit point) are derived using 
the GRB data (method of (IWang||2008l )). SNe la Union data, and BAO peak measurements, 
while dotted lines (cross denotes the best-fit point) are derived using SNela and BAO data 
only. Numerical noise is responsible for the jaggedness of some parts of contours. The best-fit 
parameters in both cases are f2m = 0.27, Q\ = 0.7 with = 321 for 307 degrees of freedom 
(dotted lines) and = 326 for 313 degrees of freedom (solid lines). Note the different axes 
scales compared to Figure 1. 




Fig. 5. — lo", 2cr, and 3a confidence level contours for the XCDM mo del. Solid li nes (circle 
denotes the best-fit point) are derived using GRB data (method of (IWandl2008l )). SNe la 
Union data, and BAO peak measurements, while dotted lines (cross denotes the best-fit 
point) are derived using only SNe la and BAO peak data only. The best-fit parameter 
values are: for solid contours (circle) - = 0.24, = —0.90 with = 327 for 313 degrees 
of freedom, and for dotted contours (cross) - f2m = 0.25, = —0.91 with = 322 for 307 
degrees of freedom. Note the different axex scales compared to Figure 2. 
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Fig. 6. — la, 2a, and 3a confidence level contours for the (/)CD M model. S olid lines (circle 
denotes best-fit point) are derived using GRB data (method of ( lWang||2008l )). SNe la Union 
data, and BAO peak measurements, while dotted lines (cross denotes best-fit point) are 
derived using SNe la and BAO data only. Numerical noise is responsible for the jaggedness 
of some parts of contours. The best-fit parameters in both cases are: f2ni = 0.24, a = 0.30 
with = 326 for 313 degrees of freedom (solid lines) and = 321 for 307 degrees of 
freedom (dotted lines). Note the different axes scales compared to Figure 3. 
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Fig. 7. — la, 2a, and 3a con fidence level con tours for the ACDM model from the GRB data, 
derived using the method of ISchaeferl (120071 ). The cross indicates best-fit parameter values 
f2m = 0.91, fl\ = 0.0 with = 77.86 for 67 degrees of freedom. The dashed line demarcates 
spatially flat ACDM models. 
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Fig. 8. — la, 2(7, and 3a con fidence level con tours for the XCDM model from the GRB data, 
derived using the method of ISchaefeii (120071 ). The best fit is achieved on the line = 0.00, 
which corresponds to the spatially flat matter-dominated Universe, = 77.8 for 67 degrees 
of freedom. The dashed line demarcates spatially fiat ACDM models. 











0.1 



0.2 



0.3 



0.4 



0.5 

n 

m 



0.6 



0.7 



0.8 



0.9 



Fig. 9. — la, 2(T, and 3a con fidence level con tours for the 0CDM model from the GRB data, 
derived using the method of ISchaeferl (120071 ). The cross indicates best-fit parameter values 
fim = 1-0, a = 4.5 with = 77.8 for 67 degrees of freedom. The a = horizontal axis 
corresponds to the spatially fiat ACDM case. 
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Fig. 10. — la, 2cr, and 3a confidence level contours for the ACDM model (dashed line de- 
marcates spatially flat models). Solid lines (ci rcle denotes the best-fit point) are derived 
using the GRB data (method of (ISchaefeii 120071 )). SNe la Union data, and BAO peak mea- 
surements, while dotted lines (cross denotes the best-fit point) are derived using SNe la 
and BAO data only. Numerical noise is responsible for the jaggedness of some parts of the 
contours. The best-fit parameters are fl^a = 0.27, = 0.7 with = 321 for 307 degrees of 
freedom (dotted lines) and fl^a = 0.28, = 0.69 with = 401 for 376 degrees of freedom 
(solid lines). Note the different axes scales compared to Figure 7. 
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Fig. 11. — la, 2a, and 3a confidence level contours for the XCDM mo del. Solid lines (circle 
denotes the best-fit point) are derived using GRB data (method of ( ISchaefen 120071 )). SNe 
la Union data, and BAO peak measurements, while dotted lines (cross denotes the best-fit 
point) are derived using SNe la and BAO peak data only. The best-fit parameter values 
are, for solid contours (circle) - = 0.26, = —0.90 with = 401 for 376 degrees of 
freedom, and for dotted contours (cross) - fl,^ = 0.25, = —0.91 with = 322 for 307 
degrees of freedom. Note the different axes scales compared to Figure 8. 
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Fig. 12. — Icr, 2o", and 3o" confidence level contours for the (pCDM model. Solid li nes (circle 
denotes best-fit point) are derived using GRB data (method of (jSchaefeii 120071 )). SNe la 
Union data, and BAO peak measurements, while dotted lines (cross denotes best-fit point) 
are derived using SNe la and BAO data only. The best-fit parameters are = 0.24, 
a = 0.30 with = 401 for 376 degrees of freedom (solid lines) and Q^n = 0.25, a = 0.30 
with = 321 for 307 de grees of freesom (dotted lines). Note different axes scales compared 
to Figure 9. 



